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\^ ' Abstract 

^^ 

We propose a mininial description of single field dark energy/modified gravity within the 

effective field theory formalism for cosmological perturbations, which encompasses most existing 

models. We start from a generic Lagrangian given as an arbitrary function of the lapse and of the 

extrinsic and intrinsic curvature tensors of the time hypersurfaccs in unitary gauge, i.e. choosing 

as time slicing the uniform scalar field hypersurfaces. Focusing on linear perturbations, we identify 

(~| ' seven Lagrangian operators that lead to equations of motion containing at most two (space or 

time) derivatives, the time-dependent coefficients of three of these operators being determined only 

by the background evolution. We then establish a dictionary that translates any existing or future 

^ I model whose Lagrangian can be written in the above form into our parametrized framework. As 

f^ ' an illustration, we show that Horndeski's — or generalized Galileon — theories can be described, up 

^T I to linear order, by only six of the seven operators mentioned above. Finally, in order to make the 

link with observations, we provide the entire set of linear perturbation equations in Newtonian 

gauge, the effective Newton constant in the quasi-static approximation and the ratio of the two 

\l ■ gravitational potentials, in terms of the time-dependent coefficients of our Lagrangian. 

o 
m 



X 






00 



Contents 

1 Introduction 2 

2 General Lagrangian in unitary gauge 5 

2.1 Background equations 6 

2.2 Perturbations in the ADM formalism 7 

2.3 The EFT language 9 

3 Generalised Galileons 11 

3.1 Geometric preliminaries 12 

3.2 Ls 13 

3.3 L4 13 

3.4 L5 14 

3.5 Summary 15 

4 Observables 15 

4.1 Perturbation equations 16 

4.2 Modification of gravity 18 

5 Conclusion 19 

6 Acknowledgments 20 
A Lagrangian dependence on ^^^R^^K^'^ 20 
B EFT parameters for generalized Galileons 22 
C Coefficients of the perturbation equations 23 



1 Introduction 

Dark energy has now become a generic name that includes a huge number of models trying to account 
for the present cosmic acceleration [1,2]. Given their proliferation, the confrontation of such models 
with present and future cosmological data would be greatly facilitated by an effective approach that 
can mediate between observational data and theory. Ideally, such a phenomenological approach 
would provide an effective parameterisation that minimizes the number of free functions and deals 
directly with the relevant low-energy degrees of freedom, which in our context are the cosmological 
perturbations (together with the background evolution) . A precise dictionary rephrasing the various 
models into this common language would then simplify the confrontation with the data and point out 
possible degeneracies between different theories. Within its unifying picture, this effective approach 
should have the extra virtue of stimulating theorists to study previously unexplored regions of the 
parameter space which could lead to interesting new models or, conversely, to better understanding 
why certain regions might be forbidden. 

A few steps in this direction have been undertaken recently. The so-called effective field theory 
(EFT) of cosmological perturbations is a powerful tool that allows to deal directly with the relevant 
low-energy degrees of freedom of the problem at hand. Such an approach was proposed and intensively 
used for inflation [3,4], in particular to characterize high-energy corrections to slow-roll models and 
to predict high-order correlation functions (see e.g. [5,6]). The EFT of inflation has now become 
a standard way of parametrising primordial non-Gaussianity and was used, for instance, in the 



interpretation of the most recent WMAP [7] and Planck [8] data. This approach has also been 
applied to dark energy, first in the minimally-coupled case [9] where it was proven a useful tool to 
study the stability in full generality and, for models with vanishing sound speed, the clustering of 
dark energy down to very nonlinear scales [10]. 

More recently, the EFT formalism has been extended to dark energy with non-minimal cou- 
plings [11, 12], providing a unifying theoretical framework for practically all single-field dark energy 
and modified gravity models.^ This approach relies on two basic steps [11]: a) assume the weak 
equivalence principle and therefore the existence of a metric g^u universally coupled to all matter 
fields (it is straightforward to relax this assumption, but at the price of complicating the formalism); 
b) write the unitary gauge action, i.e. the most general gravitational action for such a metric compat- 
ible with the (unbroken) spatial diffeomorphism invariance on hypersurfaces of constant dark energy 
field. 

In [11] it was argued that the EFT of dark energy has all the virtues advocated at the beginning 
of this section. The goal of this article is to provide a systematic procedure to translate an arbitrary 
dark energy model into the EFT language, as well as to establish a firm minimal setting of Lagrangian 
operators within this framework. In particular, here we focus our attention on the operators of the 
unitary gauge action that lead to at most two derivatives in the equations of motion for linear scalar 
perturbations.^ This minimal set of operators encompasses most of the theoretical models of dark 
energy and/or modified gravity discussed in the current literature. 

The key ingredient of our derivation is a 3 + 1 decomposition d la ADM, where time slicings 
coincide with the uniform scalar field hypersurfaces. With this time choice, the dynamics of the 
underlying degree of freedom is einbodied in the dynamics of the 3-dimensional metric. In Sec. 2 
we consider a generic Lagrangian given as an arbitrary function of the lapse N = l/\/—g^'^ and of 
the 3-dimensional metric /i^^ = g^^ + n^riy^ where n^ is the unit vector perpendicular to constant 
time hypersurfaces, more specifically of its extrinsic and intrinsic curvature tensors, respectively 
K/ = hfVpnt" and ^^^R''^, 

S = I d'x./^L{N, K^^, K^,K^\ (%, (3)i^^,(%^^ ...-t). (1) 

In our construction we include combinations of these 3-dimensional objects without taking their 
derivatives. This automatically prevents the appearance of higher (inore than two) time derivatives 
in the equations of motion. However, it is not enough to also remove higher spatial derivatives. By 
expanding this Lagrangian up to quadratic order in the cosmological perturbations and making use 
of an ADM analysis in unitary gauge (see for instance [3,18,19]) we obtain specific conditions that 
ensure the absence of higher spatial derivatives in Sec. 2.2. 

Moreover, we also show how the parameters in front of the standard EFT operators of [11] can be 
expressed in terms of the time-dependent coefficients of the expansion of (1). Since the action of most 
of the existing theoretical models can be written as eq. (1), this can be used to derive a dictionary 
between theoretical models and our EFT language. As an illustration, in Sec. 3 we explicitly derive 
this dictionary for the most general scalar field theory leading to at most second order equations 
of motion, i.e. the Horndeski theory [20] (see also [21]), recently rediscovered in the context of the 
so-called Galileon field [22,23] under the name of "generalized Galileons" [24,25]. 

Let us summarize here the main results of Sees. 2 and 3: 



^An alternative formulation of a background independent effective approach to dark energy and modified gravity 
was given in [13]. A covariant EFT of cosmological acceleration was developed in [14] for inflation and generalized to 
the case of dark energy in [15, 16] — see discussion in [11] for a comparison between the latter approach and the one 
advocated here. 

^This is sufflcient to ensure that we only have a single propagating degree of freedom. Note, however, that higher 
time derivatives do not lead to higher degrees of freedom if they can be treated perturbatively, i.e. evaluating them 
using the lower order equations of motion [17]. 



• The most general EFT action, up to quadratic order, for single-field dark energy, in the Jordan 
frame, leading to at most second-order equations of motion for linear perturbations can be 
written as 



^* fu^n ^u^ .m.oo , ^IW,. oo^2 ^m 



-f{t)R-A{t)-c{t)g^^+::^{6g^'r " '^SK6g'^ 



S = I d x^—g 



m. 



ml{t) 



(2) 



• 



where 5g^^ = g^^+1, 6K^,y = K^^—Hh^y, K = K^^ and we have assumed a flat Universe so that 
^^'R^v vanishes on the background. This action describes the propagation of one scalar degree 
of freedom with dispersion relation cj^ = c^k'^, where Cg is the sound speed of fluctuations given 
by eq. (35) with the relations (47). Stability (absence of ghosts) is ensured by the positivity 
of the time kinetic term given in (47) . The particular combination appearing in the operator 
proportional to m\ is such that it does not lead to higher-order spatial derivatives. One can 
check that also the combination 

(%^5i^;;- ^(%5i^ (3) 

does not generate higher derivatives. However, this operator is not explicitly included in eq. (2) 
because it can be reexpressed in terms of the others (see App. A). 

In the particular case where m\ = m\, the above action is equivalent to the linearized Horn- 
deski's theory/generalized Galileons. The explicit dictionary between generalized Galileons and 
this action is given in App. B. 

Expanding the Lagrangian (1) up to quadratic order we also find three operators that lead to 
higher order space — but not time — derivatives. These are 



5'h.s.d. = I d X 



.rnl{t)5K^ + ^^^^^R5K +^^^)r^ 



(4) 



When one of these operators is present in the action the dispersion relation of the propagating 
mode receives corrections at large momenta, w^ = c^A;^ -|- A;^/M^, where M is a mass scale. 
These corrections may become important in the limit of vanishing sound speed, such as in the 
model of the Ghost Condensate [26] or for deformations of this particular limit [3,9]. 

Once a Lagrangian describing matter has been included, the action (2) can be used as a benchmark 
for the study of physical signatures of dark energy/modified gravity in the linear regime. In this 
context, cosmological perturbations are usually discussed in Newtonian gauge, which is the one that 
we employ in Sec. 4. In order to do that, in Sec. 4.1 we restore the covariance via the Stueckelberg 
trick [3, 4, 26] and we vary this action with respect to all the scalar dynamical degrees of freedom. 
This allows to derive Einstein's equations and the evolution equation for the fluctuations of the scalar 
field responsible of the acceleration, recovering and generalizing previous results [9, 11]. Finally, in 
Sec. 4.2 we use these equations to derive the effective Newton constant and the ratio between the 
two gravitational potentials. 

The phenomenology of the operators appearing in action (2) was also studied in [12]. In this 
reference it was indeed mentioned that these operators are sufficient to describe linear perturbations 
of Horndeski's theories. However, no proof of this statement was given nor the particular combination 
in which such operators appear shown. 



2 General Lagrangian in unitary gauge 

In the presence of a scalar field (p with a non- vanishing timelike gradient, the so-called unitary 
gauge corresponds to a choice of time slicing where the constant time hypersurfaces are uniform 
hypersurfaces. The use of unitary gauge accomplishes two main objectives. First, as explained at 
length in Refs. [3,4,26], it makes it straightforward to write a generic Lagrangian for cosmological 
perturbations. Since the dynamics of the scalar field has been "eaten" by the metric, the most generic 
Lagrangian is simply that for the metric perturbations around a FLRW solution, compatible with the 
unbroken symmetry of 3-dimensional diffeomorphisms. 

Second, the 3+1 splitting in unitary gauge easily allows to keep the number of time derivatives 
under control, while considering higher and higher space derivatives. Therefore, the unitary gauge 
is helpful to systematically explore the space of higher spatial derivative theories by considering 
geometric invariants on the 4> = constant hypersurfaces. In practice, we will use the metric in the 
ADM form 

ds^ = -N^dt^ + hij {dx' + N'dt) {dx^ + N^dt) , (5) 

where the 3-dimensional metric hij is used to lower and raise latin indices i,j,--- = 1,2,3. Since 3- 
dimensional diffeomorphism invariance is preserved in unitary gauge, it is natural to write operators 
(with up to two spatial derivatives per field) in terms of the extrinsic and intrinsic curvatures K^i, 
and ^^'Rfiu and their possible contractions. The Lagrangian is also an explicit function of the lapse 
function A'^ in general. 

Therefore, in the following, we consider a general action of the form 

S= f d^x^/^L{N, K,S,TZ,Z;t) , (6) 

where the Lagrangian L is an arbitrary function of N and of the following four scalar quantities 
constructed by contracting the extrinsic and intrinsic curvature tensors: 

K = Kf'^, n = (3)i? = (^^i?'; , S = K^^K^"', Z = ^'^%^'-^^R^"' . (7) 

Although one should also allow, in principle, for a dependence on 3^ = ^^'R^j^yK^'^ ^ we have preferred 
not to include it explicitly in the main body, for simplicity. As shown in App. A, this extra dependence 
leads to a quadratic Lagrangian of the same form as that found later in this section, with slightly 
modified coefficients. Indeed, since y is equivalent to HTZ at linear order, the quadratic terms in the 
expansion of the Lagrangian induced by its dependence on y are analogous to those induced by its 
dependence on TZ. As for the linear term, one can use the equality 

m^^^.K^'^ = ^^'^RK + ^^""^R + boundary terms , (8) 

which is also shown in App. A. 

Moreover, one could also consider scalars that are combinations of three or more tensors, like 
K^ K^uK^x, but it is easy to show that also in this case they can be re-expressed in terms of the 
above combinations, plus corrections which are at least cubic in the perturbations. We will show 
this explicitly for the extended Galileon in the next section. Finally, one could take quadratic 
combinations of the Riemann tensor such as ^^'R^yp^^'^'R^^f" . However, in three dimensions the 
Riemann tensor can be expressed in terms of the Ricci scalar and tensor.'^ Thus, at quadratic order 



This can be done using the relation 



(3)n — (3)n u _ (3)n u _ (3)r „ + '^^R h —-^^^R(h h - h b ) CQl 



in the perturbations, the action above seems to exhaust all the possibilities compatible with our 
requirements. 

In order to explicitly write the expansion of the action (6) up to second order in the perturbations, 
it is useful to define the tensors 

5K = K- 3H, 6K^, = K^, - Hh^, , (10) 

which vanish on the background, and to use the decompositions 

5 = 3^2 + 6S, SS = 2H5K + SK'^JK"^ . (11) 

The quantities TZ and Z vanish on the background and are therefore already perturbative (-H is even 
a second order quantity). 

The expansion of the Lagrangian up to second order in the perturbations yields, after discarding 
irrelevant boundary terms, the expression 

L{N,K,S,TZ,Z)=L-T- 3HT + {f + Ln) SN + Ln STZ 

+ -Lnn Sn^ + B SK6N + C 5K5n + Lnti 5N5n + LzSZ + 0(3) , 

where we have introduced the following notations for some combinations of the partial derivatives of 
the Lagrangian (denoting Lj\[ = dL/dN, etc.), to make this expression more compact: 

T = 2HLs + Lk , 

A = iH^Lss + AHLsK + Lkk , ,^^. 

B = 2HLsN + Lkn , 

C = 2HLs-R + Lkti ■ 

The first term, L, is the homogeneous Lagrangian and all partial derivatives of L that appear in the 
above expression are evaluated on the homogeneous background, i.e. for N = 1, S = 3H^, K = 3iJ, 
TZ = Q and Z = {). Note that, in order to obtain the expression (12), we have rewritten the term 
linear in 5K as 

F5K = T{K - 3H) , (14) 

and integrated it by parts using K = V^n^, 

d^xy/^TK = - I d^xV^n^W^T =- I d*xy/^^ , (15) 

where n'^ is the unit vector orthogonal to constant time hypersurfaces and, in unitary gauge, has 
time component rfi = 1/A^. 

2.1 Background equations 

For the background we assume a fiat homogeneous FLRW metric, written in the form 

ds^ = -N'^{t)dt^ + a^{t)6ijdx'dx^ . (16) 

In this case K = 3H/N and S = 3H'^ /N'^, where H = a/ a is the Hubble rate. Note that it is crucial 
to explicitly keep the lapse function A^, because the first Friedmann equation is the constraint arising 



from the invariance under time reparametrization. Linear variation of the homogeneous action Sq 
with respect to the lapse A^ and the scale factor a yields 



6So 



dt(rx 



[L + Ln- 3HT) 6N + 3a^ [L- 3HT - J") Sa 



(17) 



where we have used \/—g = a'^N. Then the Friedmann equations are directly given by^ 

3HT -L-Ln = 0, (18) 

which depends on first time derivatives at most, and 

T + 3HT -L = 0, (19) 

which determines the dynamics of the scale factor. 

As expected, by using the above equations one can check that the first order of the total La- 
grangian C = ^ —gL vanishes. Indeed, using \J—g = yhN, where h is the determinant of the 
3-dimensional metric hij in the ADM decomposition, one easily finds 



Ci = (l - 3HF - P\ 8Vh + a^{LN + L - 3HF)5N + a^Ln^n, 
where the last term is a total derivative and can be ignored. 



(20) 



2.2 Perturbations in the ADM formalism 

In this sub-section we perform the analysis of the perturbations in unitary gauge and by using the 
ADM form of the metric, eq. (5). For the action at second order, we will only need to take into 
account the perturbations of \f—g at first order, 5^J—g = 5yh + a?5N , because the second order one 
multiplies the LHS of eq. (19). Thus, the quadratic Lagrangian for perturbations is given by 



Co =6Vh 



+ a^ 



{T + Ln)SN + Ltz 6TZ 

1 \ 1 

Ln + -Lnn I SN^ + Ln52n + -A 6K^ + B 6K5N + C 5K5n 



(21) 



1 



+Ls 5K\i 5K + Lz Sni 5ni + i^Lnn ^W + {Ln + L 



Nn, 



5N6n 



where 621^ denotes the expansion of TZ at second order in the perturbations. 

In the ADM decomposition (5) the only relevant components of the extrinsic curvature tensor 
are given by 

Kij = ^{kj - ViN, - V,Ni) , (22) 

where Vj stands for the covariant derivative associated with the 3-dimensional metric hij. For explicit 
calculations in unitary gauge we choose to describe scalar perturbations of the spatial metric in terms 

ofC[18], 



hi 



a^{t)e'^^5ij . 



Thus, the perturbations of the quantities used above are given by 

5Vh = 3a^C , SK'j = U - H5n) b) - ^6'''d^kNj) , 



and 



6n,, = -5,,d\ - didjC , 6211 = -- [{dcf - ^Cd^C] 



(23) 

(24) 
(25) 



We have not included explicitly the matter in the Friedmann equations, but it is straightforward to do so. 



Note that in this section d stands for a spatial derivative and 9^ = did^. By using the above 
expressions, the variation of £2 with respect to 5N yields the Hamiltonian constraint, which reads 



[Lnn + 2L7V + ^H {3HA + 2HLs - 2B)] 6N + 3{B- 3HA - 2HLs) C + K^n + -^)C 
-{B- 3HA - 2HLs) -^-A{Ln + Lnu - ^HC) ^ = . (26) 



By varying £2 with respect to the shift 

Ni = ditp 

one obtains the momentum constraint, which implies 



{B - 3HA - 2HLs) 5N + {A + 2Ls) 



d'^Tp 



{3A + 2Ls) C - 4C 



a^c 



(27) 



(28) 



By combining the two constraints, one can express both 6N and d'^ij) as functions of C, and its 
derivatives and then substitute in the action to write it only in terms of C, and its derivatives. In 
general, a term proportional to (i9^C)^ will remain. Here, in order to single out the lowest derivatives 
operators first, we want to find conditions under which this term disappears. If one considers the 
second order action before the substitution of the constraints, one finds the following terms 



{A + 2Ls) (aV)' + 4C aV a^C + 2 {^Lnn + 3Lz) {d^Cf 



(29) 



Taking into account the momentum constraint (28), one immediately sees that imposing the three 
conditions^ 

A + 2Ls = Q, C = 0, 4L7^7^ + 3L2 = , (30) 

implies the elimination of the term proportional to (9^C)^ ™ the final action and the absence of 
higher derivatives in the equation of motion for C,. As we will see in the next section, all generalized 
Galileon models satisfy the three conditions (30). 

When (30) are satisfied, the momentum constraint reduces to 

4L5 



6N = VQ 



V 



(31) 



B + AHLs 

By direct substitution into £2 and after an integration by parts to get rid of the term C9^C (note that 
the CC term vanishes because of the background equations of motion), we finally get the following 
Lagrangian for Q: 



A = ^ 



^aC + ^d,(,d,Q 



(diC) 



21 



with 



C 



cc 



-Lnn + Ln 



C 



9»ca.C 



L-R 



^ d , ... 
— riaM) 



3HB - QH^Ls 1 V^ + I2L5 , 

M = V{Ln + Lnu) ■ 



(32) 



(33) 



Classical and quantum stability (absence of ghosts) requires that the time kinetic energy is positive 
(see, e.g. [3,9]), 

£^^ > . (34) 

The sound speed (squared) of fiuctuations can be simply computed by taking the ratio 

2 _ -^e.C^^C 



C 



(35) 



CC 



^Note that these conditions are only sufficient. A more general analysis can be performed by explicitly requiring that 
the coefficient of (9^C)^ vanishes once the two constraints have been solved. However, this leads to a very complicated 
equation involving many of the coefficients of the quadratic Lagrangian and it is not clear whether one can find physically 
relevant solutions that evade (30). 



2.3 The EFT language 

We are now going to express the conditions on the absence of higher derivatives in terms of the 
coefficients of the action of the EFT formaUsm of Refs. [11, 12]. The action up to quadratic order in 
the perturbations can be written as 



S = d^Xy/—g 



(36) 

where R in the first term inside the bracket is the four-dimensional Ricci scalar. Note that, in order 
to make the comparison with the previous subsection simpler, we have found more convenient to 
use the 3-dimensional Ricci scalar and tensor in the quadratic terms, instead of the four-dimensional 
ones used in Ref. [11], since the link with the ADM decomposition is then transparent. 

Let us first discuss how the background equations (18) and (19) translate in this language. In 
action (36) we have used the time-time component of the inverse metric g^'^ and its perturbation 
in the expansion of quadratic and higher order operators, as it is customary in the EFT formalism. 
However, in the previous subsections it was more convenient to work directly with the lapse function 
TV, related to 5°° by 

^" = -]^ • (37) 

Only the first three terms in brackets in eq. (36) contribute to Z, Ljy and T ^ and thus to the 
background equations of motion. Using eq. (37) and employing the decomposition of the four- 
dimensional curvature scalar, 

R = (3)i? + K^^K^"' -K^ + 2V^{n^V^n>' - n^V ^n^) , (38) 

after an integration by parts in the action we can rewrite these terms as 

Lo = ^(jn + fS-fK^-2f^^-k + ^, (39) 

(we remind the reader that TZ = ^^>R). By expanding at linear order in 5N , integrating by parts the 
terms linear in X, we can match the background and linear terms of this action with the first line of 
eq. (12), which yields 



L-T - ?,HT = M^{3fH^ + 2fH + 2fH + /) + c - A 
J- + Ljv = M^ifH - 2fH - f) - 2c . 



2.... .... . . (40) 



From these two relations and using the background equations of motion (18) and (19) one finds that 
c and A are given by 

c + A = :iM^(fH^ + fH) , (41) 



K-c = M^(2fH + 3fH^ + 2fH+'f) . (42) 

This coincides with what was found in Ref. [11] in the absence of matter. 

To discuss linear perturbations we only need the second-order expansion of the action (36). By 
rewriting the first three terms as in eq. (39), expressing g^^ in terms of N and using the definitions 



(10) and (11), one immediately sees that the EFT action is of the form (6). One can thus use the 
second-order expansion of the Lagrangian (21) with the following dictionary: 

-Lnn + Ln = c + 2M| , 

A = -Mlf - Mi , 
B = JMl - fh\ , 

C = f, (43) 

Lntz = fJ-i , 
Lnn = Ai , 

which is completed with eq. (40). 

With these relations, the conditions for the absence of higher derivatives, eq. (30), can be written 
in the EFT of dark energy language. They read: 

Ml + Ml = Q, m5 = 0, 4Ai + ^A2 = 0. (44) 

These conditions are straightforward to verify. Using eqs. (24) and (27), bK"^ contains a higher 
derivative term, (S^V)^; while 5K%5K^ contains (didjip)'^. However, when the first condition in 
eq. (44) is satisfied the combination of higher derivative terms in eq. (36) gives an irrelevant boundary 
term. The second condition implies that the operator ^^'R6K, which contains d'^ipd'^C,, does not 
appear. Finally, (3)i?2 ^ iQ{d^Qf/a'^ and (^^Ri/^W^ = [5{d'^C)^ + {didjC)'^]/a'^: one can check that 
when the third condition is satisfied the sum of the two operators in eq. (36) vanishes up to a total 
derivative. 

In summary, the most general EFT Lagrangian which does not generate higher derivatives in the 
linear equations for the perturbations is 

L = ^fm - A(t) - c(t)<7°° + ^{6g'y - 4^ 5K5g'' 

where 

ml = fhl, ml = -{Ml - Mi) , ml = ^il , (46) 

as in eq. (2). Terms containing ^^'B? and ^^'Rt, ^^'R^ do not appear because they only contain higher 
spatial derivatives. By employing the dictionary (43) in eq. (33), the quadratic action for ( is given 
by eq. (32) where 



£ ■ ■ = 2 ( c + 2Ml - SH'M'^f - SHM'^f + 3/7m^ - GH^ml) V + 6(M^/ + 2mi) 



^d.Cd^C. 



Mlf-^iaM) 
a at 



(47) 
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and 

2(M2/ + 2ml) 



V 



2H{M^f + 2ml) + M^f - m^ ^48) 

M = '^{M^f + 2ml). 

The stability of a given model is then determined by the condition C^^ > and the speed of sound 
can be straightforwardly computed from eq. (35) by using the relations above. One can check that 
these results agree with those found in [11] in the limit m| = m| = 0. Note that all these relations 
hold also in the presence of matter universally coupled to the metric g^i, , except that in this case its 
homogeneous energy density and pressure pm and pm must be included in the background equations 
(see eqs. (95), (96)). 

Finally, we can also write down the independent operators that generate higher spatial derivatives. 
These are 

Lh.s.d. = - mlit) 6K' + ^ (^^R 6K + ^(3)^2 (49) 

We now turn to study a well known example of scalar-tensor theories of gravity which does not 
generate equations of motion with higher derivatives and, when restricting to linear perturbations, 
is contained in the Lagrangian (45). 

3 Generalised Galileons 

In four dimensions, the most general scalar tensor theory having field equations of second order in 
derivatives is a combination of the following generalized Galileon Lagrangians [20,24,25], 

L2=G2(</.,X), (50) 

L^=G-i{ct),X)Uct>, (51) 

L4 =Gi{(t>,X)R - 2G4x(<^,X)(n0' - </'''^>;m.) , (52) 

L5 =G^{<i),X)G^,(l)-^^^ + iG5x(<^,^) (□</>' -^U(t)(l),^,ct>'^^ + 2(/);^,</.;^-,/.'';j . (53) 

For notational convenience, in this section we mostly indicate covariant differentiation with a semi- 
colon symbol, i.e. ;^. Moreover, we have defined X = (p'^cp-fj, (note that X is sometimes defined 
differently, i.e., with a factor of —1/2). 

In order to translate the above Lagrangians into our EFT language we will proceed in two steps. 
First, we will rewrite each of these Lagrangians in terms of 3-dimensional geometrical objects {K'^ , 
^^'R'^i^, etc.) so that their unitary gauge expression becomes easily readable. The 3+1 decomposition 
that we are after loses manifest general covariance but shows straightforwardly the lack of higher time 
derivatives already at the level of the action. The second step will be to compute the corresponding 
coefficients of the operators (45) by simply inverting the dictionary that we derived in the previous 
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section — eqs (40) and (43): 
1 



J" + Ln) + HL-ji - 2L'jiH - L-ji , 

A = -L - -Ln + 2-^ + 3^-^ + '^HLn + QH^Ln + 5/7L7^ + ^Te , 

/ = 2LnM-^ , 
^1 = \ {Lnn + ^Ln+T)-]^ {HLn - 2HLn - In) , (54) 

m^ = 2Lti — 2HLsN — Lkn = 2Lti — B , 
^1 = 1 Us - 2Ln - 2H^Lss - 2HLsk - \lkk\ =\iLs- 2Ln) - \a , 

rhl = Lnu , 

where we have directly adopted the notation (46) which holds in absence of higher derivatives — more 
general relations are easily found when eq. (44) is not satisfied. 

The main result of this section is that the dynamics of linear perturbations for all generalized 
Galileons is described by (45), with the further restriction m\ = m\. This is in agreement with the 
result [24,25] that also higher space derivatives are absent from the equations of motion. This section 
is rather technical; the reader uninterested in the details of the calculations can skip the following 
subsections and go directly to Sec. 3.5 where we summarize our main results. 

3.1 Geometric preliminaries 

In order to express in unitary gauge terms of increasing complexity, it is useful to review the geo- 
metric formalism adapted to the 3 + 1 decomposition and separate the quantities into "orthogonal" 
and "parallel" to the hypersurface (j) = const. First, we define the future directed unitary vector 
orthogonal to the hypersurface. Up to a factor 7, this is proportional to the gradient of , 

The metric induced on the (j) = const, hypersurface is h^y = n^n^ + g^^^. Orthogonal to n^ are 
also various quantities that "live" on the hypersurface, in the sense that they vanish when contracted 
with n^: the extrinsic curvature and the "acceleration" vector 

K^y = /i^ n^-a, n^ = n" n^.^ . (56) 

The last two equations can be inverted by decomposing the derivative of n^ into parallel and paral- 
lel/orthogonal components, 

nu■,^Ji = K^^ - n^hu . (57) 

By means of the quantities just defined, we can decompose the second derivative of the scalar 

field as 

-1 7^ A 

4>;i^u = -7 {K^,u - n^n^ - n^n^) + — </>' XAn^n^, . (58) 

Again, this decomposition into parallel and orthogonal quantities is useful when calculating compli- 
cated products such 4>-,^u4>''^" 4>'^-.a that appear in L5, see eq. (53). 

Other relevant equations are the Gauss-Codazzi equations, relating the Ricci tensor and scalar 
intrinsic to the hypersurface, ^^'R^u and ^^'R, to the four-dimensional ones [27,28], 

^'H- = (^M-)ll + in''riPR^^,p)\\ - KK^, + K^.A'^,, (59) 

(3)i? = R + K^- K^.uK^"' - 2{Kn^ - n^);^ , (60) 

where the symbol || means projection on the hypersurface of all tensor indices, e.g. (V^)|| = h^V^. 
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3.2 Lg 

Since L2 is trivial, following [11] we start from L3 and see how to rewrite it in the EFT of dark 
energy formalism. First, it is convenient to define an auxiliary function F3{(j),X) such that 

G3 = F3 + 2XFsx . (61) 

Thus, L3 in eq. (51) can be written as 

Ls = F3ncb + 2XFsxDcl). (62) 

We integrate by parts the first term on the right-hand side and we rewrite the second term using 
Ocj) = —'y~^K + ^(l)''^X-fj,/X, which is obtained by tracing eq. (58). This yields 

L3 = -{FsxX,^ + F3^ ^■,^)<P'^ - 2Xj-^F3xK + FgxX^,/.''^ . (63) 

After noticing that the first term inside the parenthesis cancels with the last one we finally obtain 
an expression for L3 which is of the form of eq. (12), 

L3 = -2-f~^XF3xK - XFs^ . (64) 

In unitary gauge 4>(t,x) = (/>o(t), which implies, for instance, 

F3x(.^,X) ^ Fsx{Mt),-'Pl(.t)/N^) . (65) 

Using eq. (54), it is now straightforward to derive the corresponding EFT parameters in terms of 
the Lagrangian parameters evaluated on the background. They are explicitly given in App. B and 
coincide with those given in [11]. They only depend on four Lagrangian parameters, Gs,/,, G^x, 
Gsxtf) and G^xx, so that the dependence on the auxiliary function F3 disappears. As expected from 
eq. (64), /, ml and ml all vanish: in order to describe L3 we only need c. A, M2 and m^. 

3.3 U 

We now proceed with L4, defined in eq. (52). Using eq. (58) and its trace we can rewrite this as 

-1^ + ^cp-'^X.,,y - j-HK.uK^" - 2n,hn - ^(<^^'^X^)2J ^^^^ 

GiR + 2XGix{K^ - K^^K^^") + 2GixX,^{Kn^' - ii^") , 



-L4 — G4R — 2tr4Jf 



where in the second line we have used that 7 ^ = —X. Moreover, for the last term we have 
replaced 7~^(/)'^ by — n'* and used h^ = ^^h^^ X-^y. In this last term we can employ that G^xX-^ = 
d^Gi — Gi^(j)fj_ = d^G^ + 7~^G4(^n^. After an integration by parts this yields, using n^n^ = 0, 

2G4xX■f,{K7^^' - h^) = -2G4{Kn^' - n^^).^ - 27-^^4,/,^ . (67) 

The first term on the right-hand side of this expression can be rewritten by using the Gauss-Codazzi 
equation (59). Plugging all this into the second line of eq. (66) we finally obtain L4 in 3+1 decom- 
position, 

U = G4^^^R + {2XG4X - G4){K^ - K^^K^"^) - 2^~^Gi^K . (68) 

It is now lengthy but straightforward to apply our usual dictionary (54) to derive the corresponding 
EFT parameters. Their explicit expression can be found in App. B. They depend on the six La- 
grangian parameters G, G^x, G^x^)-, G^xx, G^xxcp and G^xxx- We need all the seven parameters 
of the EFT action (45) to describe L4 but the last two are equal, m| = ?fi|. 
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3.4 Ls 

This Galileon Lagrangian is more involved than the others. Let us start working on the first term 
on the right-hand side of eq. (53), G5 G^ycf)'^^ . Integrating it by parts gives 

G5 G^.<A''^" = -G5xX'>'^G^, - G^^^-^G^.n^'rf . (69) 

At this stage, as we did for L3, it is convenient to define an auxihary function F^[(j),X), such that 

G,x = ^5X + ^ , (70) 

and use this definition to integrate by parts the term proportional to G^x hi (69). In particular, 
using that 

G^xX,p = 7Vp(7~'F5) + F^^r^np , (71) 

we obtain 

G5 GpA^^^'^ = F^c^'^'^Gp, + 7"'(F50 - G5^)Gp,n>'n'' - ^F^'^'n'^Gp, . (72) 

Let us now work on the second term on the right-hand side of eq. (53). Using eq. (58) we can 
rewrite this as 



^G^xinf - 3n#;^,(/>''^- + 2(^,p,(^^^''^(t)■\) = -G^x^lC + G^xJ , (73) 



3 

where 

1C = K'^ - SKKpi.K'"' + 2Kp^K^"'K''^ , (74) 

J = -^0'^Xp(K2 _ Kp.K^^'') - 27-3 (i^n^nA* - K^.n^n'^) . (75) 

The term proportional to J on the right-hand side of eq. (73) can be integrated by parts using the 
same trick as above, which yields 

G5xJ = -F57'' Q/C + K'^^n'^nPR^^.p - Kn'^nPR^p + h'^nPR^p^ - ^^f.^k"" - K^^K^'^) . 

(76) 
For the last part of the calculation we need the (one time-)contracted Gauss-Codazzi relation, 
eq. (59), which gives 

KP'^G^, = KP''^''%u - KP'^n'^nPRp^.p + KKI, - K^, - ^RK . (77) 

Replacing (j)'P'^ with eq. (58) in eq. (72) and using this relation, the terms proportional to F^ in 
eqs. (72) and (76) combine and simplify to 

-7-^F5 (^^^Gp^KP-^ - ^JC\ . (78) 

Using this and putting together all the terms of L5 from eqs. (72), (73) and (76) we finally obtain 

i5 = - 7"'i^5 f i^'^'^^Mv - W''^^) + ^Ar^G^xIC 

V 2 ; 3 (79) 

+ -X(G50 - F^^f^R + -XG54K^ - Kp^Rp") , 

where in the last line we have used 

2Gp^n^'n'' = ^^'^R + K^ - Kp^RP" . (80) 
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Operator 


f 


A 


c 


Mi 


m| 


9 ~ 9 
mf = ml 


L2 





/ 


/ 


/ 








L3 





/ 


/ 


/ 


/ 





u 


/ 


/ 


/ 


/ 


/ 


/ 


L5 


/ 


/ 


/ 


/ 


/ 


/ 



Table 1: A list of the different contributions of the generalized Galileon Lagrangians (50)-(53) to the operators 
of (45). 

Note that the last line of (79) has the same form as the first two terms of L4 given in eq. (68): by 
using eq. (70) it can be written as 



G/^^R + (2XG4X - G4)(K' - K^.K^"') , 



(81) 



with G4 = ^X{G5^ - F50). 

In order to compute the coefficients of the various EFT operators we use the dictionary (54) . To 
treat the term K^'^^^'R^u we employ the prescription described by eq. (126) in App. A. Moreover, it 
is useful to notice that, up to quadratic order, the combination /C of the extrinsic curvature tensor 
can be replaced by an expression that depends only on S and K: 



)C = 6W 



6H^K + 3HK^ 



3HS + 0{3) 



(82) 



The EFT operator coefficients are explicitly given in App. B. One finds that they depend on the six 
Lagrangian parameters G^^, G5X, G5x<f>, G^xx, G5XX4, and G5XXX— the dependence on F5 exphc- 



itly cancels out — and, as in the case of L4, 77x4 



;^2 



m 



Thus, at linear order in the perturbations- 



quadratic in the action — L5 does not bring any new operator with respect to L4. The difference 
between L4 and -L5 appears at the cubic order in the action. 

3.5 Summary 

We have established a dictionary between the generalized Galileon theory, eqs. (50)-(53), and the 
EFT of dark energy parameters entering the action (2). Such a dictionary is explicitly given in App. B. 
As summarised in Table 1, the EFT operators and their associated time-dependent parameters that 
are needed to describe the generalized Galileons are only six: c, A and /, the three usual parameters 
already present at the background level, and Mg, m\, ni\ = m|, progressively appearing in L2, L3, 
L4 and L5, contributing only to the perturbations. As already stressed, at quadratic order in the 
perturbations, L4 contains the same number of independent operators as L5 — in particular, only 



the combination m\ 



ml 



appears in the action. The case 772.4 7^ ^4 encompasses the generalized 



Galileons: it does not contain higher derivatives and yet does not belong to the generalized Galileon 
class. When m^ 7^ ml, higher derivatives are expected to appear beyond linear order. However, the 
effect of these higher derivatives can be ignored as long as perturbations remain small and linear 
theory is a good approximation. 

4 Observables 

Observables describing large scale structures are computed in the framework of linear cosmological 
perturbation theory. In this section we first derive the perturbation equations describing the dynamics 
of dark energy and modified gravity. We include a matter sector describing cosmological species such 
as cold dark matter, baryons, photons and neutrinos — by adding the matter Lagrangian JCm{9iiu,i^m) 
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to eq. (2), so that the final Jordan frame action in unitary gauge reads 
S = d!^X\/—q 

J L ^ 

ml{t) 



^fm-m-c{t)g''+ ^{Sg'r - ^SKSg^' 



- mlit) {6K' - 6K,M<n + ^^'^Rdg'^ (83) 

-fhli 



i{t)5K' + ^(')R6K + M(3)^2 ^ c^(^g^^,^l.^ 



2 2 

We then discuss the modifications of gravity expected in hnear theory. We wih use Newtonian gauge, 
which is often used in cosmology, especially to describe cosmological perturbations for modified 
gravity. Extension to other gauges or to so-called "gauge invariant" formalisms is straightforward. 

4.1 Perturbation equations 

We will first restore the general covariance of the action above and write it in a generic coordinate 
system. In order to do that we need to reintroduce the scalar fluctuation tt via the Stueckelberg 
trick [3,4,26]. Under the time coordinate change t — t- t + Tr{t,x), the four-Ricci scalar R remains 
invariant, while functions of time such as / and the 3-dimensional quantities change as^ ^ 

f^f + h + ^fn^ (84) 

gOO ^ ^00 ^ 2^0^^ ^ g^-^d^TTd^TT , (85) 

6Kij — ^ 6Kij — H-rrhij — didjir , (86) 

6K ^5K - SHtt - ^d^TT , (87) 

^^^Rij -^ ^^^Rij + H{didj7r + 5^jd'^7r) , (88) 

i^)R^(^)R+A^Hd\. (89) 

In the new coordinates we consider a linearly perturbed FLRW metric with only scalar fluctua- 
tions, 

ds'^ = -(1 + 2^)dt'^ + 2dia dtdx' + a^{t) [(1 - 2^)5^^ + 2xij] dx'dx^ , (90) 

where Xij is traceless and given in terms of the scalar perturbation /3, Xij = {didj — ■^6ijd'^)f5. The 
extrinsic curvature and the 3-dimensional Ricci tensor of the new equal-time hypersurfaces thus read 



K,j = e~''{H-^)h,,+Xij-did, 



a 



' ^ ' (91) 

^^^Rij = didj^ + 6,jd^^ + 2dkd(,xl) - ^''x^J ■ 

We also decompose the matter stress-energy tensor at linear order as 

T''o = -{pm + 6pm), (92) 

T\ = {Pm + Pm)^^V = -o'T^ , (93) 

T) = {pm + Sp„^)S} + (d'd, - ^5pA a , (94) 

®With an abuse of notation, here we denote the extrinsic curvature on hypersurfaces of constant time with Kij even 
when we are not in unitary gauge. The reader must be aware that Kij is not the same geometrical object before and 
after the Stueckelberg trick. The same also holds for '^^^Rij. In particular, after the Stueckelberg trick Kij and ''^'-Rij 
are respectively given by eq. (91). 

^The operator rh\ is also considered in Ref. [12]. However, in vl of this reference, the variation (89) of "'_R under the 
Stueckelberg trick has been overlooked and the error propagates into the Einstein equations and the various observables. 
With the authors of [12] there is now agreement on this issue [29]. 
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where p^ and Pm. are respectively the background energy density and pressure and Spm and Sp^, their 
perturbations, v is the 3- velocity potential and a the scalar component of the anisotropic stress. The 
background equations derived from the action (83) are [11] 



A - c = Ml{2fH + 3//72 + 2fH + /) + p^ . 



(95) 
(96) 



Using these expressions and the transformations (84)-(89) allows us to rewrite (83) as an action 
for the scalar fluctuations $, a, ^, /3 and vr. We can vary eq. (83) expanded at second order and then 
fix the Newtonian gauge by setting a = = /? in the equations derived. This yields five equations, 







1 <55(2) 



1 55(2) 



A^^ + A^i> + A^^ + A^^ + \{Af^ + ^^V) + 5T\ , 



(97) 



a=0=l3 



l—g ba 

1 55(2) 



= k' 



a=0=l3 



B^^ + B^^f + S^vr + B^ir + -^(5$?^^ + ^[2)71 



ik'ST^i , (98) 



-55^ 



1 55(2) 



C$$ + C|,6 + C^j,^ + C^^ + C^vr + Cf,-fr + C^ii- 



a=0=/3 



k' 



k^ 



a"* 



■<5/3 



1 55(2) 



+ ^(Cf $ + CS?)^ + d'K + Cf ^) + ^(4^)^ + d'K) - 5T\ , (99) 



a=0=/3 



i^^fc^") k'k, (42)$ + dJ)^ + ^^2)^ ^ ^(2)^ ^ ^(^2)^ 



U' 



+ \{Df^> + D^^\)]-5T 



(100) 



-fif^vr 



£;$$ + E^^ + E^-^ + E^^ + Ei^-^ + E^TT + Ei,-k + E^^n 



a=0=l3 



.(2) 



(4), 






(101) 



The coefficients Aa, Ba,Ca, Da and E'a of these equations are detailed in App. C. 

The explicit expressions of the above coefficients given in the appendix contain also higher- 
derivative terms — those proportional to ml, fh^ and A. To compare with the usual Einstein equa- 
tions, here we rewrite these equations by replacing the components of the stress-energy tensor Tt 
with their expressions given in eqs. (92)-(94). For simplicity, we set ml = m^ = A = 0. We obtain: 



• 00-component {6S/6^ = 0): 

M2 



2 



\a 



kl 

2' 



2/ ( ^^ + ZH^ + 3//2$ ) + / ( " vr + 3//2vr - 2,H{^ - ir) - 3(^ + H^) ) + ?,Hf'K 



(c + A)7r + (2c + 4M| + ^Hml) {^ - ir) + {ml - 4.Hml) 

k^ 
4^m|(^ + Htt) = 5pm ■ 



U2 



• Oi-component {5S/5a = 0): 

Ml [{Hf - /)7r + / ($ - vi-) + 2f{H^ + "J)] - 2ctt - m| {^ - n) + 4ml{H^ + 'J + Htt) 

= - {Pm. + Pm) V . 



(102) 



(103) 
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ij-trace component (53/6"^ = 0): 



M^2f 

+ / 



3 a 



2 ($ - ^) + (3i/2 + 2i7)$ + H{^ + 3^) + ^ 



2 k' 



■K + 2iJ$ + 2H{^ --k)- {3H^ + 2H)tt + 2^ + 6 - if 



3a2 

+ / [-2H7r + 2($ - 7f )] - /(^V J> + (A - c)7r + 2c($ - ir) 
4A;2 



(104) 



3a^ 



^2 [m^($ - vf) + (Fm^ + {ml)') vr + ??t,^7t] 
- AmlH + (mlY + 3Hml ($ - tt) + 4(ii'm4)'7r - m|(6 - vf) 
+ 4i7 [(m^)" + 3Hml] $ + 4(7n^)"^ + 4mlH{3HTr + 6 + 3^) + 4m^§ = 5pm ■ 
• ij-traceless component (63/ 6/3 = 0): 

M^ /($-^) + /7r +2[mln + mlHTT + {mlyTT]+2ml{<^-ir)=a . (105) 

By combining eqs. (97) and (98) we obtain the relativistic generalization of the Poisson equation, 

1.2 

F$$ + F^^ + F^TT + F^^ + — (i^i^^^ + 42) vr) = 5pm - 3H{pm + P^m)v = PmAm , (106) 

which can be also written as: 



Generalized Poisson equation: 



k^ 



{2fMt + 4mi)^ - {fMt 



m- 



+ AHml - AHmhir 



+ (GM^H'^f - 6Hc -C-A + 3mlH)7r 



- (2c + 4M^)n + {3MtHf + 2c + 4M|)$ - SM^f'^ + 3?n^(^ + H<^) = pmA 



4.2 Modification of gravity 



(107) 



In order to derive the effective Newton constant, Gqs, we consider the quasi static approximation, 
i.e. we neglect the time derivatives in the equations of motion and we neglect the anisotropic stress, 
(T = in eq. (100). This is a good approximation for scales much smaller than the sound horizon 
scale, i.e. for k ^> aH/cg. For models with small or vanishing sound speed (see e.g. [10]) or on scales 
longer than the sound horizon, a consistent treatment which takes into account the time derivatives 
should be undertaken. 

In the quasi-static limit, GcS is defined by 



k' 



^ = 4TTGcs{t,k)pmAr 



(108) 



Following [30], in order to write the Poisson equation in this form we can use eqs. (100), (101) and 
(106). For Cs ~ 0(1), we can neglect D^, Dqr, D^, E^, Ei^/, F^ and F^ from these equations and the 
effective Newton constant is thus given by 



47rG, 



eff 



-[M-'[ 
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M = 




l(2) 



(4), 



Ef+E'i\k/af 



.(2) 



(4)/ 



E^{k/ay^+Ei^^ +Ei^\k/af 



F,: 



(2) 






(2) 



(109) 
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We can write it in a slightly more compact form as 



o_2(fc/a) ^ + Oo + 62(K/a) 



where 



0-2 — D^ Ej^ , 

a2 = i^r 4^) - 4^)i^r - 4^^4^ + ^?^4) , 

a4 = -4^)4'^ + 4'^4), (111) 

6_2 = 4'^^-4^ ' 

62 = -4^4^4) - 4^)4V(f) + 4)4^)4^ + 4^4^4) • 

Another quantity often used to parameterize deviations from General Relativity is the ratio 
between the gravitational potentials 7 = ^/<^, which is given by 

[com(.M)]32 
^=[com(Al)]3,' ^''^^ 



where com(A^) denotes the comatrix of M. This reads 

_ c_2(A:/a)~2 + co + C2(A;/a)2 



a-2{k/a) 2 + ao + a2{k/aY + a4(/c/a)'^ ' 
with 



(113) 



c_2 = -4^^- ' (114) 

CO = 4^4^ - 4^4^ , (115) 

02 = 2^(^)4^-4)4'^. (116) 

The expressions for Geg and 7, eqs. (110) and (113), generalize those given for instance in [30] 
in absence of higher derivative operators, in which case 02 = 04 = 62 = C2 = 0. When also 
o_2 = 6-2 = c_2 = we recover the results of [11]. It was stressed in [31] that the numerator of G^s 
equals the denominator of 7. Here this is confirmed also in the presence of higher derivatives.^ 

5 Conclusion 

In this paper we lay down the basic building blocks for a systematic phenomenological study of dark 
energy and its cosmological perturbations. Following [11], our basic assumptions are that a) dark 
energy/modified gravity brings in at most one scalar propagating degree of freedom and that h) the 
weak equivalence principle is satisfied — there exists a metric tensor universally coupled to matter. We 
use the effective field theory formalism developed for inflation in [3,4], that is based on an expansion 
in number of perturbations rather than in number of fields. Indeed, expanding the action in number 
of fields [15, 16] becomes unpractical every time that the background field configuration undergoes 
a large excursion. On the opposite, the main advantage of the present (non-covariant) approach is 
that an expansion in number of perturbations can always be consistently truncated at the desired 



^It simply follows from [A4 ^]i3 = (det A^) ^[com(7\/()]3i 
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order of approximation, in virtue of the empirical fact that perturbations are small on the largest 
scales. Moreover, our formalism is "ready to go", in the sense that there is no need of solving for the 
background equations first. Apart from the three operators /, c and A responsible for the background 
evolution [11], every new operator is at least quadratic in the perturbations: it does not affect the 
background and its dynamical effects can be studied independently. 

In particular, we consider only operators that are at most quadratic in the number of perturba- 
tions — those needed for the linearized equations of motion — and we single out a set of seven operators 
that bring up to two derivatives in the equations of motion. To achieve this result, in Sec. 2 we provide 
a systematic treatment of any Lagrangian that can be written in ADM form as a general function 
of extrinsic and intrinsic 3-dimensional curvature tensors and of the lapse function. This is already 
enough to avoid higher time derivatives in the equations of motion. Then, in Sec. 2.3 we identify 
specific combinations of the EFT operators that are required to avoid higher-order spatial derivatives. 
Some operators can be re-expressed into other ones, thus simplifying the EFT Lagrangians up to 
quadratic order. 

The entire Horndeski, or "generalized Galileon", theory can be written in this formalism (Sec. 3): 
a relevant amount of work has gone into re-expressing all the generalized Galileon Lagrangians in 
their ADM form and obtaining their EFT formulation. At linear order, Horndeski theories can 
be described by a total of six operators: only three quadratic operators in addition to those — /, 
c and A — accounting for the background (see eq. (2) with ml = ml). This seems a substantial 
simplification if compared to the full covariant treatment and well represents the power of the non- 
covariant EFT approach. The two Galileon Lagrangians L4 and L5, despite their scaring looks (52)- 
(53), are affordable at linear order in the perturbations with the addition of just one operator with 
respect to those needed for L3. At linear order, there exists another operator beyond the Horndeski 
theory (for ml 7^ rhf) that still gives equations of motion limited to second order in time and space 
derivatives. Finally, we analyze also some higher spatial derivative operators, those in eq. (4). 

The time dependent coefficients of our seven plus three operators described by actions (2) and 
(4) remain to be constrained or measured by observations. Indeed, in Sec. 4.1 we provide the set 
of linear perturbation equations in Newtonian gauge by varying these actions with respect to scalar 
metric and field fluctuations in a generic gauge. As an illustration, using these equations we com- 
pute the effective Newton constant in the quasi-static approximation and the ratio between the two 
gravitational potentials (Sec. 4.2). The computation of these " observables" should be considered as a 
first step towards a more general and systematic study of the impact of dark energy on cosmological 
perturbations in order to fully exploit future observational data. 
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A Lagrangian dependence on ^^^Rn^K^^ 



In this appendix we show how to treat a dependence on 

3^ = ^^^.K'^" (117) 

in the unitary gauge Lagrangian. 
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Let us first show the relation 



I d^x^^\{t)'^^'^R^uK^"' = I (fxy/^ 



^i^)RK + ^^^)R 



2N 



or, equivalently, 



j4 / 

a x\j—g 



\{tf'^R^,K^- - ^(3)^ ^ _ M(3)^ 



2N 







(118) 



(119) 



up to some irrelevant boundary terms. Since K = V^n'^, the last two terms in the above integral 
can be simplified via an integration by parts, so that the expression reduces to 



Jd'xV^xit) (^(3)i?^,K/^- + !^v^(%) 



(120) 



It remains to show that this can be written as the integral of a total derivative. 

Using the explicit expressions for the extrinsic curvature in the ADM decomposition, eq. (22), 
and n'^ = —Ng^^, the above expression can be rewritten as 



fd^xVhX{t) - [h'^y^hJ^'^R.j + (3)i?^ _ V'N^^^^ - -N'Vi^^^R 



(121) 



where Vj is the covariant derivative with respect to the three- metric hij. The second term can be 
integrated by parts and then vanishes when combined with the last term, as a consequence of the 
Bianchi identity V^^'^'Gij = 0. Finally, the term in parenthesis can be rewritten as 

h^^W^hki '-^^Rij + '■^^R = h'^y^hki ^^'^Rij + h'^^^^Rij + h'^'-^^Ry = h'^ '■^^Rij . (122) 

and it is known that the last expression can be reexpressed as the divergence of a three-vector, i.e. 
fiij ^^'Rj_j = ViJ^ (the very same property is used to derive Einstein's equations from the Einstein- 
Hilbert action^). We have thus proved Eq. (118). 

Let us now assume that the Lagrangian L introduced in Eq. (6) also contains an explicit depen- 
dence on y. By noting that y is already a perturbative quantity, i.e. vanishes in the background, 
and can be decomposed as 

y = Hn+ (3) R^,^ SK^" , (123) 

where the first term on the right hand side is a first (and higher) order quantity while the second term 
is only second order, one immediately finds that the expansion of the Lagrangian, up to quadratic 
order, will yield the following extra terms with respect to the expression (12) obtained in the main 
body: 

L{N, S, K, 7^, y, Z) D Lyy + (Lj^ybN + LKybK + LsybS + Lny^Tl) HbTZ + ^^Lyyii^bV? . (124) 

The first term can be expressed in terms of TZ and K by using eq. (118) with A = Ly. Expanding 
up to second order then yields 

Lyy = -ihy ^ 3HLy) 6n + - (Ly5K - Ly6N] 6TZ + 0(3) + boundary terms , (125) 



^See for instance Eq. (7.5.14) of Ref. 
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so that the expansion of the full Lagrangian now reads 

L{N,S,K,TZ,y,Z) =L-P-2,HF + Ln5N + ]- (2Ln + Ly + 3HLy) 571 



1 



+ Ls SKJ^SK;: + ( 2H^Lss + 2HLsk + ^Lkk) 5K^ + ^Lnn6N^ 



{Lnn + H^Lyy + 2HLyn) 6n^ + {2HLsn + Lkn) SK6N (126) 

+ (2HLsn + Lku + HLKy + 2H^Lsy + ]^Ly\ 6K6n 
+ (lnti + HLNy - ^Lyj SNdTZ + 0(3) . 

In summary, an explicit dependence of the action on 3^ can easily be included in our treatment, via 
the following substitutions in Eq. (12), 

1 ■ 3 

L-R.^ Lti + -Ly + ^HLy , 

Lnn -^ Lnn + H Lyy + 2HLyTi , 



Lntz —^ Lnti + HLNy — -^Ly 



1 , (127) 

1 



C — 7- C + HLxy + 2H Lsy + —Ly . 

B EFT parameters for generalized Galileons 

Here we explicitly give the EFT of dark energy parameters in terms of the Lagrangian (12), for the 
generalized Galileon Lagrangians eqs. (51)-(53). All quantities in the expressions below are calcu- 
lated on the background. 



/ = 0, (128) 

A = ^^{4> + 3H^)G3x, (129) 

c = <i)\-4> + 3H^)Gsx + 'j>^G3^, (130) 

Ml = y('A + ^H<p)G3x - 3H,p'^Gs,xx - yGg.x^ , (131) 

m^ = 2^^G3x , ml = ml = 0. (132) 

M2/ = 2G4 , (133) 

A = ^> + SHXG^x - im'^Gix4>'^ + QHG4x^(p^ + l2H'^G4xx<j)^ , (134) 

c = -\^T + 3HXG4X - m'^G^x^'^ + 6HG4x^<j)^ + l2H^G4xx<P^ , (135) 

Ml = \^- ^HXG^x + GHG^x^^^ + ISH^G^xx^^ - GHG^xx^^"" - UH^G^xxx^^ , (136) 



m 



3 — ovr^ .. Qun ..I'i 



2XG^x - mG^xr + 4G4X0'/'' + IQHG^xxr , (137) 

ml = ml = 2Gix^'^ , (138) 
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with 



T = 2MtHf + Mtf + F = 2XGix - SHG^x^, 






Mtj = -G5^cl)' + 2G5 



X(t 



f5XX<; 



(139) 

(140) 
(141) 

(142) 

M| = i> - ^M^Hf - ^H^G^x^P^ + GH^G.x^^'' + GH^G^xx^'' - SH^G^xx^^' - 2H^G^xxx4>' 

(143) 



c = -^T+^M^Hf - 3H^G,^<P^ - 3H^G5x4>^ + SH^G^x^^^ + 2H^G5 -^^ 



li- 



r2 tt2 



2/~<„ 12 ^u'in a3 I Qu2/-Y___ 14 I ouS/ 



A = -^ + 3iV4^F^/ + -MtHf - 6if^G5^0^ - IH-'G^xr + SH'G^x^^ + ^H'G^xx 



^3 7\#2 ^ 



m^ = Mtf - AHG^^cj)' - QH'G^xr + ^HG^x^p^ + ^H^G^xx^ , (144) 

ml = ml = G5^</>2 + HG^x'P^ - G5X<P^ , (145) 

with 



T = 2MtHf + Mtf + F = 2MtfH + M^f - 2HG^sr - 2H'G^x 



2 ou2r'_. 13 



(146) 



C Coefficients of the perturbation equations 



In this appendix we define the coefficients appearing in eqs. (97)^(101). For convenience we have 

(147) 



used M4 defined as 



Mi = 2mi + 3mi 



Variation v^ith respect to $: 

A^ = 2c + 4M| - 6H \fHMt + M^/ - m| + HM^ 
A^ = -3 [2H [fMl + Ml) + Mlf - mi] , 

A^ = SH'^M^f + 3mlH - c - A - QMIHH + SM^/ 
A^ = -2c - 4M| - 3H{ml - M^f) , 

4'^ = -2/m2 + GHm^ - Ami , 

^(,2) =M^f_^3^ 2HMl - AHml + GH^fh^ . 

Variation w^ith respect to a: 

B^ = -ml + 2H [fMl + Ml) + M^f , 

B^ = 2 {fMl + Ml) , 

B^ = -2c + 2MlH + M^{Hf - f) , 



B^ = ml- Mlf 



B 



(2) 



3 
* - -2^5 , 
^(.2) = -2{ml + Hm^) 



(148) 

(149) 

(150) 
(151) 
(152) 
(153) 

(154) 
(155) 
(156) 
(157) 
(158) 
(159) 
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Variation with respect to ^: 



C$ 



II 



2c + 2(3i/2 + H)M\ 

+2fMl{2,H^ + 2H) - (mi)- + H (-3mi + AMlf + 2(M|)-) + 2Mlf 
C^ = -3mi + QH [fMl + M|) + SM^/ , 
C4, = 6 {2,fHMl + 3HMi + M^/ + {Uiy 
C^ = 6 (/M2 + M|) , 



Ctt 



c - A - 2H{Mly - 2{H + 3HH)Mi + M^{2fH + /^^^ + SiJ^ + 2ii'/) 



C^ = 3 (-2c + 3Hml - 2HM^f + 2M|iJ + (mi)" - 2M^f 

Ci = 3(mi - Mlf) , 
4'^ = - (2/M2 - 6i/m5 + 4ml) , 
4^ = 2/Af 2 - 6Hm5 - 6^5 , 
C^2) = _2 (^M^/ + (M|)- + HMi + 3/^-2^5 + SHrh^^ , 



Cl 



(2) 



-2(M| - 2ml + 3Hfh5 



a 



(4) 



16A 



Ci^^ = -2m5 + 16if A . 
• Variation ^vith respect to (3: 



D 
D 



D 



(2) 
(2) 
(2) 



Mlf + 2m2 _ 3^5ii- , 
-3m5 , 



L'l^) = M^/ + 2m2i/ + 2(m2)- - 3iJm5 , 



D 
D 



(2) 

TV 

(4) 



2m4 — 2m4 , 
-8A, 



^(4) = mg - 8i7A 



(160) 
(161) 
(162) 
(163) 
(164) 

(165) 
(166) 
(167) 
(168) 
(169) 

(170) 
(171) 
(172) 

(173) 
(174) 
(175) 
(176) 
(177) 
(178) 
(179) 
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Variation with respect to n: 



E^ = 6cH + c + H^{9ml - GM^f) + 3(2mi - M^f)H -A + 3H 4M| - 2M|i7 + (ml) 



E^ = 2c + 4M| + 3i7(mi - M^f) 



E-^ 



E^ — 3 



QcH + c + k- 3M^f{2H'^ + H 

2c + 3Hml - AHM^f - 2M|if + (m; 



3) 



E, 
E. 



^ - 3{mi - Mtf) , 



GMJH'^ - 3{mlyH + 6Hc - 9HHml + c - 3M^Hf - GH'^M^f - 3mlH + A 

Ei = -2 [3H (c + 2M|) + c + 2(M|)"] , 
E^ = -2 (c + 2M|) , 



n(2) 



ml + H (-2M| - Gii'ms + 4m^) - M^/ 



K 



(2) 



2Hml + M^V - SmgiJ + 2( 



7714 1 



^f ^ = 2M| + 6iJ?fi5 - 4m^ , 



i?p) 



2-2 



™2^ 



2c - 4M|if + 4m|iJ + {ml)' + 4i/^m| + iJm^ - Um^HH + AH {mi)' 



e: 



(4) 



^ - -2m5 + 16HX , 
£;(.^) = -2(ml + 2Hfh5) + lei/^A . 

• Generalized Poisson equation: 



F$ = 3M^Hf + 2c + 4M| + 3Hml , 

F^ = -3Mlf + 3ml , 

F^ = GM^H'^f - 6Hc - c - A + Sm^H 

F^ = -(2c + 4M|), 

.7(2) 



-i2fM^ + Aml) 

t2 _3 



F^2) ^ y^2 _ ^^3 ^ 4^^2 _ 4^^2 ^ 



+ 4(M|) 
(180) 
(181) 
(182) 

(183) 

(184) 

(185) 

(186) 
(187) 

(188) 

(189) 

(190) 

(191) 

(192) 
(193) 

(194) 
(195) 
(196) 
(197) 
(198) 
(199) 
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